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Statistical isotropy of primordial perturbations is a common assumption in cosmology, but it is an 
assumption that should be tested. To this end, we develop cosmic microwave background statistics 
for a primordial power spectrum that depends on the direction, as well as the magnitude, of the 
Fourier wavevector. We first consider a simple estimator that searches in a model-independent 
way for anisotropy in the square of the temperature (and/or polarization) fluctuation. We then 
construct the minimum- variance estimators for the coefficients of a spherical-harmonic expansion of 
the direction-dependence of the primordial power spectrum. To illustrate, we apply these statistics 
to an inflation model with a quadrupole dependence of the primordial power spectrum on direction 
and find that a power quadrupole as small as 2.0% can be detected with the Planck satellite. 
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I. INTRODUCTION 

It is well known that the homogeneity and isotropy 
of the Universe are only approximate. There are de- 
partures from homogeneity and isotropy that are now 
well-quantified by measurements of the cosmic microwave 
background (CMB) and galaxy surveys. In current 
cosmological theory, the notions of homogeneity and 
isotropy have been superseded by the notions of statis- 
tical homogeneity and isotropy. The density of matter 
may differ from one point in the Universe to another, 
but the distribution of matter is described as a realiza- 
tion of a random field with a variance that is everywhere 
the same and the same in every direction. This is gener- 
ally the prediction of structure- formation models, and in 
particular, of inflationary models. 

Still, statistical isotropy and homogeneity are assump- 
tions that can be tested quantitatively, and the precision 
with which they can be tested is improving rapidly with 
the still-accumulating wealth of cosmological data. Pre- 
liminary (and controversial) indications for a preferred 
direction in the CMB [J H, H, 0] have recently motivated 
the study of departures from statistical isotropy. Subse- 
quent theoretical work has shown that although statisti- 
cal isotropy is a generic prediction of inflation, inflation 
models can in fact be constructed to violate statistical 
isotropy @, 1, 0, & i M HI M 01 M ■ Dark-energy 
models might also accommodate departures from statis- 
tical isotropy [HI, [l(| . These models provide useful straw 
men against which the success of the standard inflation- 
ary predictions of statistical isotropy can be quantified. 

The growing interest in such models motivates us to 
study generalized tests for statistical isotropy. In a statis- 
tically isotropic Universe, the primordial distribution of 
matter is a realization of a random field in which Fourier 
modes of the density field have variances, a power spec- 
trum -P(k), that depend only on the magnitude k of the 
wavevector k. If we drop the assumption of statistical 
isotropy, the power spectrum will depend on the direc- 
tion k as well. If S(k) is the Fourier amplitude of the 



fractional density perturbation, then the power spectrum 
is defined by 



{S(k)S*(k')) =S D (k-k')P(k), 



(1) 



where the angle brackets denote an average over all re- 
alizations of the random field, and Sd is a Dirac delta 
function; note that we are still preserving the assump- 
tion that different Fourier modes are uncorrelated. The 
most general power spectrum can then be written, 



P(k) = A(k) 



J29LM(k)Y LM (k) 

LM 



(2) 



where YlmW (with L > 2) are spherical harmonics, and 
<7lm(&) quantify the departure from statistical isotropy 
as a function of wavenumber k. Since the density field is 
real, Fourier modes for k are related to those of — k, in 
such a way that the multipole moment L must be even. In 
the limit <7lm(&) — > 0, we recover the usual statistically 
isotropic theory with power spectrum A(k). The imple- 
mentation, Eq. @ , of power anisotropy is motivated in 
part by the inflationary model of Ref. [10J , which predicts 

52JVf(fc) 7^ 0. 

Here we consider several CMB tests for statistical 
isotropy. The first, which we refer to as "power mul- 
tipole moments," is a simple and intuitive estimator that 
involves measurement of the multipole moments of the 
square of the temperature/polarization fields. 1 As an 
example, we apply this statistic to an inflationary model 
[Tfl | that predicts a quadrupole in the matter power spec- 
trum. 



1 There has already been some evidence for a dipole in the CMB 
power 0, 3 that is analogous to the higher multipole mo- 
ments that we are considering here, but which cannot be due 
to anisotropy in the primordial power spectrum because it has 
L = 1. There have also been searches Q| for anisotropy along the 
lines considered here, and Ref. [17, 18] discusses similar statistics. 
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Although power multipole moments provide a nice 
model-independent test for departures from statistical 
isotropy, more sensitive probes can be developed if the 
particular form of the departure is specified. To illus- 
trate, we thus construct the minimum-variance estima- 
tors for the anisotropy coefficients gLM{k) under the 
assumption that they are constants. The naive power 
multipole moments, although intuitively simple, co-add 
a number of modes with equal weight. The minimum- 
variance estimator co-adds modes with weights that de- 
pend on their signal to noise, so that (as the name sug- 
gests) the variance of the estimator is minimized. We 
show that this statistic provides a far stronger probe for 
the gLMS. 

The plan of our paper is as follows: Section HT1 reviews 
some CMB basics. Section ITO1 calculates the correlations 
of CMB spherical-harmonic coefficients if there are de- 
partures from statistical isotropy. As we discuss there, 
the power spectrum C;, which describes the two-point 
CMB statistics if there is statistical isotropy, is general- 
ized to a set of moments Df^ 1 if statistical isotropy is 



broken. In Section IIV1 we introduce and calculate the 
power multipole moments and calculate the standard er- 
rors with which these moments can be recovered. We 
apply this statistic to a quadrupole in the matter power 
spectrum, calculating the sensitivities of several CMB 
experiments to such a quadrupole. Section [V] discusses 
minimum-variance estimators for the quantities Du, M 
that parametrize the departures gLAi(k) from statistical 
isotropy. We then construct from these the minimum- 
variance estimators for the quadrupole moments of the 
primordial power spectrum, calculate their variance, and 
evaluate their sensitivity to departures from statistical 
isotropy. We make some concluding remarks in Sec- 
tion IVII Throughout the main body of the paper, we 
discuss statistics for only a temperature map, in order 
to make the presentation clear. An Appendix general- 
izes to include the full temperature-polarization infor- 
mation. Our numerical results are for a full temperature- 
polarization map, as well as for temperature or polariza- 
tion alone. 



dent for different I and m: (az m a^, m ,) = Ci6u>6 mm > . 2 The 
set of C;s is the CMB temperature power spectrum. We 
will see that when statistical isotropy is violated, there 
are correlations induced between a; m s for different I and 
m [10j. If there is statistical isotropy, the two-point au- 
tocorrelation function is 



C(n,n') = (T(n)T(n')) 



(4) 



i.e., the correlation function depends only on the sepa- 
ration between the two points. If statistical isotropy is 
violated, this is not necessarily true. 



III. OFF-DIAGONAL CORRELATIONS FOR 
ANISOTROPIC POWER 

Consider a primordial matter power spectrum P(k) 
given by Eq. ([2]). We expand T(n) in k space in the 
form, 

^(n) = Jd 3 kj2(-i) l (u+m$-n)m®m, (5) 



where &i(k) is the contribution to the Ith temperature 
moment from wavevector k. With these conventions, 
Oz(fc) is real. With our expression, Eq. we can write 
the covariance matrix as 



'm') - $U'5mm'Cl + ^ tfJti'm' D U^ • ( 6 ) 



Here, the set of C;S, given by 



LM 



Ci = (4tt) 2 / dkk 2 A(k)[ei(k)} 2 , 
Jo 



(J) 



is the usual CMB power spectrum for the case of statis- 
tical isotropy. Departures from statistical isotropy intro- 
duce the second term, where 



II. PRELIMINARIES 

A CMB experiment provides the temperature T(n) as 
a function of position n on the sky. The map T(n) can 
be expanded in terms of spherical harmonics Yj m (n), 



ah, 



i- J dnY lra (n)T(n). 



(3) 



DfiF = (4ir) 2 (-i) 1 - 1 ' / dfcfc 2 A(% LM (fc)eKfc)e^(fc), 

Jo 

(8) 

and 



3^f, w , = / dk^(k)^ ro ,(k)y LM (k) 

' (r<L \ 1 / 2 r <LM 



(9) 



The ai m s are Gaussian random variables, and if there 
is statistical isotropy, then they are statistically indepen- 



2 Strictly speaking, it is not the a; m s that are statistically inde- 
pendent, but rather their real and imaginary parts. 
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where CkN, m , are Clebsch-Gordan coefficients, and 



L _ (2I + 1)(2I' + 1) 



4tt(2L + 1) 



(Qoi'o) 



(10) 



Throughout, we use upper-case indices LM for power 
anisotropics, and lower-case indices Im for tempera- 
ture/polarization anisotropies. For L even, £i m l'm' are 
nonvanishing only for I — I' even, and so the Dfp are real. 
Eqs. dU) and j9|) agree with similar results in Ref. [l8j |. 
and they recover the results of Ref. [l(| for L = 2. 

If primordial perturbations are statistically isotropic 
and Gaussian, then the statistics of the CMB tempera- 
ture map are specified fully by the power spectrum, the 
set of C;S. If primordial perturbations have a departure 
from statistical isotropy that can be written in terms of 
spherical harmonics YtAf(k), then the two-point statis- 
tics are described additionally by the set of multipole 
moments Df(, . These quantities are thus the general- 
ization of the Cis if there is statistical anisotropy. 



IV. POWER MULTIPOLE MOMENTS 

A. Theoretical predictions 

It is natural to expect that a spherical-harmonic pat- 
tern of anisotropy in the matter power spectrum mani- 
fests itself in a similar pattern in the CMB power. It is 
thus natural to consider a set of "power multipole mo- 
ments," 



bhM — 



i J dnY* LM {h) <T 2 )(n), 



(11) 



where (T 2 ) (fi) = C(n, n) is the expectation value of the 
square of the temperature at position fi in the sky; it 
is the autocorrelation function at zero lag. With this 
statistic, we simply look for anisotropies in the power. 
These statistics have several advantages. In addition 
to having a form familiar from similar statistics [e.g., 
Eq. ([3|)] for temperature fluctuations, they have simple 
analytic expressions in terms of -P(k). There are also 
(as we show below), relatively simple expressions for the 
cosmic-variance- and instrumental-noise-induced errors 
in the measurement of these statistics. 

The variance (T 2 ) (fi) as a function of position fi is 
given by 



(T 2 ) (fi) 



T 2 



E 

Iml' m' 



\ai m al ml )Yi m {h)Y{; ml {n). (12) 



We put this into Eq. (fTTj) and use 



^Iml.—m' ^ImV 



L'M' 



mm' 



Sll'Smm 1 , 



to obtain (for L > 2) 



w 



(13) 



(14) 



B. Statistical noise 

We now calculate the standard error, due to cosmic 
variance and instrumental noise, with which the power 
multipole moments can be measured. To do so, we con- 
sider a full-sky map T map (fi) of the temperature in -/V p i x 
equal-area pixels. The temperature in each pixel re- 
ceives contributions from signal and from noise. Thus, 
in pixel i, T ma P = T(n;) + where T(hi) is the tem- 
perature measured in pixel i, which will be the signal 
temperature smoothed by a Gaussian beam of full-width 
half maximum (fwhm) 0f w hm 7 plus a noise Tf. We as- 
sume that the noise is isotropic and that the noises in 
different pixels are uncorrelated with variance a\: i.e., 
^T"T") = a\bij. The power spectrum for the map is 
thus C, map = \Wi\ 2 C t + Cf, where Cf = {A-k /N v ^)a 2 T is 
the noise power spectrum, and Wi is a window function 
that takes into account the effects of beam smearing; for 
a Gaussian beam of fwhm 0f w h m , it is Wi = exp(— l 2 a\ /2) 
with a h = 6> fwhm /V8 In 2 = 0.00742(6» fwhm /l°). 

Since the instrumental noise is isotropic by assump- 
tion, we get an unbiased estimator for &lm (for L > 2) 
from 

KTi =^sj dnY£ M (h) [T map (n)f . (15) 

Cosmic variance and instrumental noise induce a variance 
in the &lms, which wc define as 



1 U — ( bLM^LM 



(16) 



where we have assumed the null hypothesis, qlm = 0. 
For this null hypothesis of a statistically isotropic Gaus- 
sian random map, 



■-LM 



[ dndn'C map (n,n')C map (n,fi') 



xY LM (n)Y£ M (n') 

2 ~y ' qL pmapj-imap 
W 



(17) 



where C map (fii, n.2) is the two-point correlation function 
for the map, obtained from the expression, Eq. for 
the correlation function by replacing C\ by C ; map , and we 
have used J2mm'(^tmi'm') 2 = !• Note that the absence of 
any M dependence of S^A^ i s as we expected. Moreover, 
it follows from Eq. (|13[) that the estimators for the dif- 
ferent bi m s are uncorrelated: ^ijf^i'M' / ^ Sll'Smw- 

Given a power spectrum of the form Eq. |2]) , specified 
by the functions gLAi(k), predictions for the b L f% can be 
evaluated with Eq. (fl"4|) replacing Dhf^ in that equation 



byD w 



r\LM 

W 



WiWf and evaluating the Df/ 1 with 



Eq. ijHJl. The 6™^ p can then be measured using Eq. 
with variances given by Eq. (jTTJ) . 
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C. A worked example 

As a simple example, suppose the gLM{k) are con- 
stants, independent of k. We can then take gLM outside 
the integral in Eq. ((8]). An estimator for gLM is then 

g LM = bZV{bZVg L M). Defining F w = D w M /g LM 
for this case, the variance with which each gLM can be 
measured is then 



[E«< G^,F lv WiWv 



9lm 



(18) 



Moreover, the measured gLM are statistically indepen- 
dent as a consequence of the statistical independence of 
the b L M- 

To illustrate, we apply this result to an inflationary 
model [To| that has a power spectrum with a quadrupole 
dependence on the angle. 3 We use the &i(k) calculated 
by CMBFAST 19] to obtain F w . We assume only scalar 
perturbations and the current best-fit cosmological pa- 
rameters. 

The numerical results are given in Table 1, but before 
reviewing them, we provide some very rough estimates 
to get some feel for the numbers. To do so, ignore instru- 
mental noise and suppose that Wi — 1 for all I < l mayi - 
For L = 2, F w ^ only for V = I or V = l± 2. Moreover, 
for these combinations of IV and for Z 3> 2, we approx- 
imate the numerical results (which we use for the nu- 
merical results in the Table) for Fu> as i*y+2 — —0.5 Cj. 



IS 



Also, (C?° o y ~ (5/8)*- 1 for / » 1, and (C™ 1±2)0 ) 
1.5 times as large. Eq. (fT5|) can then be approximated 



rrn, - 2567rE^(C*z) 2 ]/E^C* i ] 2 . If the power spec- 



7 2 

giM 



trum has the form Ci cx I 2 (a very rough approximation 
to the temperature power spectrum for I < 1000), then 



2 

^min 2 and L a v 



a 



128ttC 



hi(^maxAmin)] ■ For example, using 
= 1000 yields a g2M ~ 1.23. 
Of course, there is nothing about the derivation of 
Eq. I|18p that is specific to a temperature map, and this 
result can be applied equally well, e.g., to the E-mode 
polarization. If we approximate the polarization power 
spectrum by Ci ~ const, then we find & g2M — 512nl~l^, 

0r Vg%M ~ 5 X 10_2 for ^max ^ 1000. 

We now return to the numerical results for a g2M listed 
in Table 1 for the Wilkinson Anisotropy Probe (WMAP) 
(20| . which has now collected three years of data, the 
Planck satellite [U, to be launched in 2008, and EPIC 
[22[, a satellite mission currently under study. The pa- 
rameters assumed for each model are listed, as well as re- 
sults obtained using Eq. (|14|) assuming only TT is used or 
EE only. The Appendix generalizes Eq. (fl4|) to the case 
where the full temperature-polarization is used (includ- 
ing the TE correlation), and we present numerical results 



3 Note that our 320 i s (2/3)v/47r/5g„, where g„ is the coefficient 
in Ref. [lfj of (k • z) 2 if the preferred direction is taken to be z. 



for this case in the Table as well. We also list results, la- 
beled "CVO" (cosmic variance only), for a hypothetical 
experiment that has perfect angular resolution and no 
instrumental noise. These numbers are for hypothetical 
full-sky experiments, but a realistic experiment will likely 
only be able to use ~ 65% of the sky for cosmology. If so, 
then each estimate for a g2M must be increased by a fac- 
tor (0.65) -1 / 2 , about 25%. We also note that the theory 
cannot specify the direction e of the quadrupole, and so 
a search for a quadrupole would require evaluation of all 
five <72AfS. A "3cr" detection would thus require that the 
sum of the squares of the 32MS need to exceed (3<7 S2JU ) 2 , 
which is independent of M. 

The order of magnitude that we would expect for o g2M 

is ~ -^Vpj,^ 2 , where -/V p j x ~ i^ ax is the number of resolu- 
tion elements on the sky, comparable to the precision 
with which one can measure the variance (the monopole) 
of the temperature-fluctuation amplitude. The numeri- 
cal results listed in Table 1 for the error to giM obtained 
from the power quadrupole moment b%M are not quite as 

— 1/2 

good as this N pix ' expectation. The origin of this dis- 
crepancy can be traced to two sources. First of all, the 
two-dimensional CMB signal is degraded from the three- 
dimensional power spectrum; a Fourier mode in the z 
direction gives rise to some temperature fluctuation near 
the north pole, and not just at the equator. This is man- 
ifest in the large coefficients (e.g., the factor of 512-71") in 
our analytic estimates. 

However, another reason that the estimator &2M does 
not provide a sensitive probe of a quadrupole departure 
from statistical isotropy is that it is not an optimal esti- 
mator for g2M- This estimator sums the "signals" -D^V , 
but it does not weight these signals properly. This can 
be seen by noting that for a C; oc l~ 2 power spectrum, 
for example, the error obtained from Eq. (|18p can be re- 
duced by applying a low-pass filter: i.e., by increasing the 
minimum values of IV in the sums. (A simple calculation 
shows that with the properly chosen lower-^ limit, a g2M 
can be reduced by a factor of 30.) If the precision of the 
result is improved by removing data, then something is 
sub-optimal. 



V. THE MINIMUM- VARIANCE ESTIMATOR 

A. The estimator and its variance 

The bLM estimator is a simple and intuitive quantity 
that can be measured to test for statistical isotropy in 
a model-independent way. However, if one has a spe- 
cific theory, defined by the functions gLM(k) or some 
quantities that parametrize the gLM(k), then there will 
be estimators that can be constructed to measure opti- 
mally those parameters. For example, if the 3la/s are 
all constants, then one can measure them better than 
the numerical results for the power multipole moments 
bLM would suggest. Below, we will derive the minimum- 
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f), V. 


pmm /rnrps 


pmm mm 
°S2M 


pmm f f f n 


mv f'T"T'1 
°92M I 1 L ) 


rr mv (F,F,1 


rr mv ('total') 
°92M l LOLcll J 


WMAP 


30.0 


42.6 


21' 


1.3 


11 


1.2 


0.024 


2.4 


0.024 


Planck 


13.1 


26.8 


5' 


1.6 


0.16 


0.16 


0.0052 


0.033 


0.0050 


EPIC 


0.021 


0.068 


52' 


1.2 


0.55 


0.42 


0.016 


0.019 


0.011 


Cosmic variance 











1.8 


0.014 


0.014 









TABLE I: The standard error <7g 2M to the amplitude of a quadrupole anisotropy in the matter power spectrum for different 
experiments. The instrumental temperature and polarization noises and beam width are listed for each experiment. We show 
results for the power multipole moments (pmm)for TT only, EE only, and the full result. We also show in the last three columns 
CT 92A/ fr° m the minimum- variance estimator for each experiment, for TT only, EE only, and the full result. 



variance estimator for qlm- 

Before moving on, it is instructive and will be useful 
below to re-derive the variance to 6 £ JJ P . We return to 
Eq. (fT4"|) and note that can be written as a sum 



M .map 

over D u , v 



We then return to Eq. ^ 



to derive the minimum- variance estimator for £)^, ' map . 



Gmap 
iven a map a 

for D 



LM ,map 



lm , each mm' pair provides an estimator 
through 



with variance 



map map,* s~t £ X 
LM, map _ a l m a l'm' ~ U l°U'°mm' 



/j^LAf,map\ 2 \ _ (1 + $XV $mrn' ) C™ aP C™' 
\ IV ,mm' J 



(cLM > 
\SXmVm' i 



(19) 



(20) 



these quantities, and we have their variances. For a 
theory with constant c/lm's, each DM* provides an es- 
timator through gLM,w = / Fw . We then sum 
these, inversely weighted by their variance to obtain the 
minimum- variance estimator, 



T. V >iFwW l WvD 



LM .map 
W 



Qlm 



( p; LM, map \ ' 

\ ) 



\LM .map 



obtained from the entire map. The variance a], of 
this estimator is then obtained by summing the inverse 
variances of all the estimators. Again, the sums are over 
/' > I, but the factor (1 + 5 lv ) in Eq. (J22J) allows us to 
write the sum over all 



The estimators for different mm 1 pairs are uncorrelated 
(if we use the real and imaginary parts of the a; m s), so the 
estimators can be summed over all mm' pairs, inversely 
weighted by the variance, to obtain a minimum-variance 
estimator. If I = I' , we sum only over to' > to to avoid 
double-counting pairs. However, the factor (1 + Sw S mm >) 
then weights the m = m' modes twice as much, if I = I' ', 
and thus allows us to re-write the sum over all to and 
to'. The result for the estimator can thus be written, for 
both I = I' and I ^ V , as 



piLAf,map 

XV ~ 



map map,vIM 
Hm "'I'm' Simi'm' 

U XV 



(21) 



We recognize these to be the bipolar-spherical-harmonic 
coefficients of Refs. pjj, with a slightly different weight. 
The variance of this estimator is then 



D 



LM, map 
XV 



(i + s, 



;;/)c7 map c map 



(22) 



The variance, Eq. (|18[) . with which each 6/ m can be mea- 
sured simply follows by summing the variances of each 
term in Eq. (fT4]) . 

Now, to construct the minimum- variance estimator, 
we simply note that the statistically-independent quan- 
tities predicted by the theory are the Dili's, the gen- 
eralizations of the C/'s for a theory without statisti- 
cal isotropy. We have constructed above estimators for 



1 



QLM 



w 



{FwWiWy) 2 



(24) 



B. Illustration: The Power Quadrupole 

To illustrate, we now evaluate this expression for L = 
2. Again, in this case, the only IV combinations that 
contribute are V = I and I' = I ± 2. We assume 1,1' 3> 1, 
approximate -Fu+2 — — 0.5 C;, as above, and evaluate 
C;gf,Q as in Section HV CI We can then write, 



1 



32 jU 



: 0.035 



E 



(C, map ) 2 



(25) 



which we can further approximate as 0.017 / max , where 
l max is the multipole moment at which Cf ~ Ci(Wi) 2 . 
The end result is then a g2M ~ 7.6/Z max , quite close to 
what we would have expected by simply counting the 
number A p j x ~ l max of usable pixels. For the WMAP and 
Planck temperature maps, Z max is roughly 650 and 2000, 
respectively, implying a g2M ~ 1.2 x 10~ 2 and 3.8 x 10~ 3 , 
respectively, implying very significant improvements in 
the sensitivity over the power multipole moments. 

Table 1 lists the exact numerical results, obtained 
by evaluating Eq. (f!?4")) exactly, for both TT only and 
EE only. Again, the Appendix generalizes Eq. (|24p for 
the full temperature-polarization map, including the TE 



6 



cross-correlation, and numerical results for this case are 
also included. The Table shows that by weighting the 
modes correctly, we get an improvement of a factor of ~ 2 
for WMAP and Planck EE and more than an order-of- 
magnitude improvement for WMAP and Planck TT; this 
is in accord with our arguments that the signal-to-noise in 
the TT power multipole moments was particularly poorly 
chosen. Although EPIC will have vastly improved instru- 
mental sensitivity, with its modest angular resolution, it 
is not particularly well suited to search for departures 
from statistical isotropy. Again, the minimum-variance 
numbers in the Table must be increased by about 25% 
to account for partial-sky coverage. And again, since the 
preferred direction is not known a priori, the sum of the 
squares of the f/2AfS must exceed (3<r g2M ) 2 to claim a "3a" 
detection of a departure of statistical isotropy. 



VI. CONCLUDING REMARKS 

We have considered CMB tests for the statistical 
isotropy of the primordial power spectrum. The power 
spectrum of Eq. (J2J) is the most general power spectrum 
if the assumption of statistical isotropy is dropped. In 
the more general case, the CMB power spectrum C; is 
generalized to a set of moments D^,, which are closely 
analogous to the bipolar-spherical-harmonic coefficients 
of Refs. 



17] . The power multipole moments &lm pro- 



vide simple and intuitive statistics that can be used to 
search in a model-independent way for departures from 
statistical isotropy. If, however, a particular model is 
introduced by specifying a particular parametrization of 
the functions <7lm(&), then minimum- variance statistics 
can be introduced to improve the precision with which 
these parameters can be constrained. For example, we 
constructed explicitly the minimum-variance estimators 
for the coefficients gLM for the case in which they are 
^-independent. We applied these results to a model in 
which there is a quadrupole in the primordial power 
spectrum, and the results are shown in Table 1. We 
see that the best probe of a primordial quadrupole mo- 
ment will come from Planck TT, for which we antici- 
pate (Jg 2M = 0.0052. Multiplying this by 1.25 to account 
for a 65% sky coverage, and then by the factor of 3 re- 
quired for a "3fi" detection, we find that the smallest 
quadrupole amplitude that will be detectable by Planck 
will be around 2.0%. 

To reduce clutter in the equations and to keep the main 
line of reasoning clear, we have derived equations in the 
main body of the paper for the case where either the 
temperature or the polarization is used, but not both. 
The Appendix generalizes the analysis to allow the use of 
the full temperature-polarization information, including 
the TE cross-correlation. 

What about other probes? Consider, for example, the 
Sloan Digital Sky Survey (23|. The volume and galaxy 
density of the main galaxy survey allows measurement, 
roughly speaking, of the amplitudes of N mo d 



dependent Fourier modes of the density field, in the lin- 
ear regime, and these measurements are cosmic-variance 
limited. Measurement of the quadrupole of the power 
spectrum can then simply be done by comparing the am- 
plitudes of Fourier modes in different directions. The 
standard error to the power multipole moments will thus 
be o~ gLM ~ <J 2 /N mo dcs ~ 10 -2 , comparable in order 
of magnitude to what can be achieved with the CMB. 
Of course, a realistic search will be hampered by the ir- 
regular volume of the survey, redshift-space distortions, 
and anisotropics (line-of-sight-versus-angular) inherent 
to the measurement technique. But then again, there 
will be degradations (foregrounds, sky cuts, etc.) to the 
idealized CMB measurements we have considered. Of 
course, if gLhiik) varies with k, then the constraints pro- 
vided by the CMB and galaxy surveys will be comple- 
mentary, to the extent that the wavenumbers k probed 
by the CMB and galaxy surveys differ. Looking forward, 
there is ultimately the possibility of accessing with 21-cm 
fluctuations approximately 10 15 modes of the primordial 
density field [2J|, allowing values as small as gLM ~ 10 -7 
to be probed, but this is in the very far future. 
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APPENDIX A: GENERALIZATION TO A 
TEMPERATURE-POLARIZATION MAP 

For most experiments, the sensitivity to departures 
from statistical isotropy will come primarily from cither 
the temperature or the polarization. Considering both 
in tandem will provide some improvement in the result, 
but given the temperature-polarization cross-correlation, 
this improvement will be weaker than what would be ob- 
tained by simply adding the two results in quadrature. 

Still, to be complete, we include expressions for theory 
and estimators for a combined temperature-polarization 
map. Assuming only primordial density perturbations 
contribute to the temperature-polarization map, a map of 
the sky will now provide the E- mode polarization E(n), 
constructed in the usual fashion (25l . [26[ from the mea- 
sured Stokes parameters Q(n) and U(h), in addition to 
the temperature T(n). The map can be written in terms 
10 5 in- of spherical-harmonic coefficients a^L, for X = {T,E}, 
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and Eq. ([6]) is generalized to 



LM.XX' 



„X X',*\ _ r r r XX' , V^tUI n L, 
a lm a i> m ' / — °U'°mm'^l t / Mml'm ,1J lV 

LM 

(Al) 

The Cf x s and ' s are obtained as in Eqs. (|7|) 
and © by replacing the 0;(fc)0/'(fc) factors in the inte- 
grands of those equations by Qf{k)Qf (k), where these 
are obtained from Eq. ([5|) by replacing T(n) by X(h). 

Note that for TE and I ^ V, Z^, M ' XX ' ^ D^ M,XX '. This 
will affect the equations below for the minimum- variance 
estimator. 

We now have a set of three power multipole mo- 
ments b^ M , obtained from Eq. (fTTj) by replacing (T 2 ) 
by (XX'), which is itself obtained from Eq. (|12p by us- 

m S lii*) f° r ^ ne expectation value therein. The 

expression for the b£ M is the same as Eq. (TT4")) using 

D LM,XX' , 

The power- multipole- moment estimators & XX /' map 
are as in Eq. (TJ5D with [T map (n)] 2 replaced by 
[X ma P(n)x' ma P(n)]. Things get a bit trickier, though, 
when we calculate the variances, as the estimators for 
different XX' will now be correlated, although still un- 
corrected for different LM. The variance in Eq. p^|) is 
now promoted to a 3 x 3 matrix S^^, for {A, A'} = 
{TT, EE, TE}. For {A A'} = XX' = {TT, EE}, 
-.if, is given by Eq. |T7| with c™ p C™ ap replaced by 
CA,ma Pc ,A',ma P _ F or the diagonal TE-TE term, 

~LM _ \ ^ f~ih r^TT.map^fEE.map . ^.TE.map^TE.mapn 
"TE,TE — 2-^i H' 1 I I' + I I' J> 

IV 

(A2) 

and for the off-diagonal XX- XX' terms, 

r^LM 2 \ , QL £-tXX,map£-<XX',map (A3) 

L ^XX,XX' / j 11' I [' \ ) 



Eq. (|18p for the standard error with which a constant 
9lm can be recovered with the power multipole moments 
is then replaced by [H, HE [27| 



1 



AA' 



<% A '* 
UU LM 



dg 



LM 



(A4) 



This is the equation used to obtain the "total" results 
listed in Table 1 for the power multipole moment. 

The minimum- variance estimator for §lm and its 
variance are similarly generalized. The estimators 
D,^, M ' A ' map are still uncorrelated for different IV pairs and 
different LM, but they are now correlated for different 



A. The main subtlety is that since 

^iv ' + D vi > 
we must be careful to keep track of all TE modes for 

I =/= V . This will require that we split the sum in the 

generalization of Eq. into two sums: the first over 

1 = 1', and the second over V > I. (Actually, the sum 

can in fact be written over all IV , but at the cost of much 

uglier algebraic expressions.) 

For V = I, there are now three (TT, EE, and TE) 

estimators to replace that in Eq. (p~9|) . and for V > I, 

there are now four (TT, EE, TE, and ET) estimators 

to replace that in Eq. (Til?]) . For all IV, the estimators 

are as in Eq. (fT!))) , replacing each a)" ap and C ; map by the 

appropriate a X ^ map and C XX ,map j respectively. The esti- 
mator for each IV , obtained after summing over all mm' , 
is the same as in Eq. (|21[) . For I = V , the variances 

^n' M ' A maP ) \ are now promoted to a 3 x 3 covari- 

ance matrix, and for V > I, they are promoted to a 4 x 4 
covariance matrix. In both cases, the covariance matrix 
can be written as 



(( 



Ml' 
C AA' 



I f^LM, A, map f^LM, A', map\ / * r\ 



For any IV pair, the diagonal entries, for A = {TT, EE}, 

are Cf A = cf' map Cy' map , and the TT-EE off-diagonal 

entry is C^t ee = C'; TE ' map C i T E ' map . For I = V, the diag- 
onal TE-TE' entry is 



--TE,TE 



^tTT^map^tEE.map 



/2. (A6) 



For /' > I, we have C# E , TE = C^C/P, C^ T ET = 
Cj T Cf E , and C l 4' E ET = C™C™. For any IV, we have 



nil' _ pTT^TE j nli 

L "TT,TE ~~ °i Z' 111111 ^EE.TE 



HI' 



C EE C ; TE . For /' > /, 



we also have C 



ii' 

TT,ET 



= CPCr and C 



u' 

EE.ET 



The generalization of Eq. ([24| is then 



^E G "E^ A '(^) 4 [(c")" 



A A' 



Y, G u'Y, F u' F iviWiWi > 



i>>i 



AA' 



c"' 



AA 
-1" 



AA' 

(A7) 



where the matrix inversion is in the 3x3 AA' space in 
the first sum and in the 4x4 AA' space in the second 
sum. We use Eq. (|A7|) to evaluate the standard errors for 
the "total" minimum-variance estimators listed in Table 
1. 
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